Specular Andreev reflection in graphene 
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By combining the Dirac equation of relativistic quantum mechanics with the Bogohubov-De 
Gennes equation of superconductivity we investigate the electron- hole conversion at a normal-metal- 
superconductor interface in graphene. We find that the Andreev reflection of Dirac fermions has 
several unusual features: 1) The electron and hole occupy different valleys of the band structure; 

2) At normal incidence the electron-hole conversion happens with unit efficiency in spite of the 
large mismatch in Fermi wave lengths at the two sides of the interface; and, most fundamentally: 

3) Away from normal incidence the reflection angle may be the same as the angle of incidence 
(retro-reflection) or it may be inverted (specular reflection). Specular Andreev reflection dominates 
in weakly doped graphene, when the Fermi wave length in the normal region is large compared to 
the superconducting coherence length. We flnd that the transition from retro-refiection to specular 
reflection with decreasing doping is associated with an inversion of the voltage dependence of the 
subgap conductance of the interface. 



PACS numbers: 74.45.-l-c, 73.23.-b, 74.50.-|-r, 74.78.Na 



The interface between a superconductor and a metal 
may reflect a negatively charged electron incident from 
the metal side as a positively charged hole, while the 
missing charge of 2e enters the superconductor as an 
electron pair. This electron-hole conversion, known as 
Andreev reflection [l|, is the process that determines 
the conductance of the interface at voltages below the 
superconducting gap — because it is the mechanism 
that converts a dissipative normal current into a dissi- 
pationless supercurrent. By studying the reflection of 
relativistic electrons at a superconductor, we predict an 
unusual electron-hole conversion in graphene (a single 
atomic layer of carbon, with a relativistic energy spec- 
trum 0, ^)- While usually the hole is reflected back 
along the path of the incident electron (retro-reflection), 
the Andreev reflection is specular in undoped graphene 
(see Fig. ^. We calculate that the difference has a 
clear experimental signature: The subgap conductance 
increases with voltage from 4/3 to twice the ballistic 
value in the case of retro-reflection, but it drops from 
twice to 4/3 the ballistic value in the case of specular 
reflection. 

The practical signiflcance of this investigation rests on 
the expectation that high-quality contacts between a su- 
perconductor and graphene can be realized. This ex- 
pectation is supported by the experience with carbon 
nanotubes (rolled up sheets of graphene), which have 
been contacted succesfully by superconducting electrodes 
0, la, la. ■ The one-dimensional nature of transport in 
nanotubes explains why the possibility of specular An- 
dreev reflection was not noted in that context, since it 
is an essentially two-dimensional effect. From a more 
fundamental perspective, the unusual Andreev reflection 
in graphene teaches us something new about the inter- 
play of superconductivity and relativistic dynamics — 
something which was not known from earlier studies of 
relativistic effects in heavy-element superconductors Q- 

We consider a sheet of graphene in the x — y plane. 
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FIG. 1: The top two panels show reflection processes that oc- 
cur in a normal metal: specular reflection at the interface with 
an insulator and Andreev retro-reflection at the interface with 
a superconductor. Arrows indicate the direction of the veloc- 
ity and solid or dashed lines distinguish whether the particle 
is a negatively charged electron (e) or a positively charged 
hole (h) . In graphene there is a third possibility, indicated in 
the bottom panel. Like in the usual Andreev reflection, the 
electron is converted into a hole, but the refiection angle is 
inverted. 



A superconducting electrode covers the region a; < 
(region S), while the region a: > (region N) is in the 
normal (non-superconducting state). Electron and hole 
excitations are described by the Bogoliubov-De Gennes 
equation Q, 



H - Ef 

A* Ef 



A 



= e 



(1) 



with u and v the electron and hole wave functions, 
e > the excitation energy (relative to the Fermi en- 



ergy Ep), H the single-particle Hamiltonian, and T the 
time-reversal operator. The pair potential A(r) couples 
time-reversed electron and hole states. 

For X > the pair potential vanishes identically, dis- 
regarding any intrinsic superconductivity of graphene. 
For x < the superconducting electrode on top of the 
graphene layer will induce a nonzero pair potential A(a;) 
via the proximity effect (similarly to what happens in a 
planar junction between a two-dimensional electron gas 
and a superconductor [13 )• The bulk value Aqc**^ (with 
4> the superconducting phase) is reached at a distance 
from the interface which becomes negligibly small if the 
Fermi wave length A^ in region S is much smaller than 
the value Xp in region N. We therefore adopt the step- 
function model 



A(r) 



Aoe'-^ if a; < 0, 
if a; > 0. 



(2) 



We assume that the electrostatic potential U in regions 
N and S may be adjusted independently by a gate voltage 
or by doping. Since the zero of potential is arbitrary, we 
may take 
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FIG. 2: Excitation spectrum in graphene, calculated from 
Eq. JHJ with A = = [/ for two values of the Fermi energy 
Ep = Tivkp- Red lines indicate electron excitations (filled 
states above the Fermi level, from one valley), while blue lines 
indicate hole excitations (empty states below the Fermi level, 
from the other valley). Solid and dotted lines distinguish the 
conduction and valence bands, respectively. The electron-hole 
conversion upon reflection at a superconductor is indicated by 
the arrows, for the case of normal incidence {k = kx, ky — 
0). Specular Andreev reflection (right panel) happens if an 
electron in the conduction band is converted into a hole in the 
valence band. In the usual case (left panel) , electron and hole 
both lie in the conduction band. In each case, the electron- 
hole conversion happens with unit probability {\rA\ = 1) at 
normal incidence, in spite of the large wave length mismatch 
between the normal and superconducting regions. 



Uir) 



-Uo if X < 0, 
if X > 0. 



(3) 



For Uq large positive, and Ep > 0, the Fermi wave vector 
kp = 2tt/\'p — {Ep + Uo)/hv in S is large compared to 
the value kp = lix jXp ~ Ep jhv in N (with v the energy- 
independent velocity in graphene). 

The single-particle Hamiltonian in graphene is the two- 
dimensional Dirac Hamiltonian llj. 



H = 











H± = -iTiv{(7xdx ± cfydy) + U, 



(4) 
(5) 



acting on a four-dimensional spinor 

(*A+,*B+,*A-,*B-)- The indices A,B label 

the two sublattices of the honeycomb lattice of carbon 
atoms, while the indices ± label the two valleys of the 
band structure. (There is an additional spin degree of 
freedom, which plays no role here.) The 2x2 Pauli 
matrices (7^ act on the sublattice index. 

The time-reversal operator interchanges the valleys 



m, 



T 



a, 

7, 



C^T- 



(6) 



with C the operator of complex conjugation. In the ab- 
sence of a magnetic field, the Hamiltonian is time- reversal 
invariant, THT^^ — H . Substitution into Eq. (^) results 
in two decoupled sets of four equations each, of the form 



H±-Ep A 

A* Ep - H± 



= e 



(7) 



Because of the valley degeneracy it suffices to consider 
one of these two sets, leading to a four-dimensional 
Dirac-Bogoliubov-De Gennes (DBdG) equation. For 
definiteness we consider the set with H+. The 
two-dimensional electron spinor then has components 
(mi,U2) = {^ A+i^ B+)i while the hole spinor v = Tu 



has components (vi,W2) = (^* 



-m 



Electron exci- 



tations in one valley are therefore coupled by the super- 
conductor to hole excitations in the other valley. 

A plane wave (u,w) x exp{ikxX + ikyij) is an eigenstate 
of the DBdG equation in a uniform system at energy 



A\^ + {Ep-U±hv\k\y 



(8) 



with |fc| = (fc^ -t- fc^)^/^. The two branches of the exci- 
tation spectrum originate from the conduction band and 
the valence band. The dispersion relation (jH)) is shown in 
Fig. 12 for the normal region (where A = = U). In the 
superconducting region there is a gap in the spectrum 
of magnitude |A| = Aq. The mean- field requirement of 
superconductivity is that Aq <?C Ep + Uo, or equivalently, 
that the superconducting coherence length ^ = Tiv/Aq is 
large compared to the wave length A^ in the supercon- 
ducting region. The relative magnitude of ^ and the wave 
length Xp in the normal region is not constrained, and 
we will compare the two regimes Xp <C ^o f^nd Xp ^ ^q. 
Simple inspection of the excitation spectrum shows the 
essential physical difference between these two regimes. 
Since ky and e are conserved upon reflection at the in- 
terface X = 0, a general scattering state for x > is 
a superposition of the four fc^^-values that solve Eq. © 
at given ky and e. The derivative h~ de/dkx is the ex- 
pectation value Vx of the velocity in the x-direction, so 



the reflected state contains only the two A:2;-values having 
a positive slope. One of these two allowed A:a;-values is 
an electron excitation {v — 0), the other a hole excita- 
tion {u = 0). As illustrated in Fig. [3 the reflected hole 
may be either an empty state in the conduction band (for 
e < Ef) or an empty state in the valence band (e > Ep)- 
A conduction-band hole moves opposite to its wave vec- 
tor, so Vy changes sign as well as Vx (retro-reflection). 
A valence-band hole, in contrast, moves in the same di- 
rection as its wave vector, so Vy remains unchanged and 
only Vx changes sign (specular reflection). For e ^ Aq 
the retro- reflection dominates if Ep ^ Aq, while specu- 
lar reflection dominates if Ep ^ Aq. 

To calculate the probability of the electron- hole conver- 
sion, we match a superposition of states with allowed k^- 
values in N and S, demanding continuity at a; = 0. The 
calculation is described in the Appendix. We give the re- 
sults for A^ <C Xp,£,, in the two regimes Ep 3> Aq, e and 
Ep ^ Ao,e. The amplitude ta for Andreev reflection 
(from electron to hole) is 



M(e, a) 



rA{e,a) 



e *■* cos a 

(e/Ao)cosa-f C 



, if ^i. > £, (9) 



e ^"^ cos a 
e/Ao + Ccosa 



, a Ef<^ e, (10) 



while the amplitude r for normal reflection (from electron 
to electron) is 



r{e,a) 



"Csina 



(e/Ao)cosa + (' 



if Ep > e, (11) 



— (e/Ao)sina 

r{e,a - .V °' , if Ep « e. 

e/Ao + Ccosa 



(12) 



Here a is the angle of incidence (as indicated in Fig. ^J 
and C = (eVAg - 1)^/2 if e > Aq, C = «(1 - e^/A^y/^ if 
e < Ao. Notice that the two regimes of large and small 
Ep are related by the substitution e/Ag ^^ C- 



One readily verifies that |r| 



rA 



1 if e < Aq, as 



it should be since transmission into the superconductor 
is forbidden below the gap. At normal incidence (a = 0) 
we find |r^p = 1 for e < Aq, so the electron- hole con- 
version happens with unit probability. This is entirely 
different from usual normal-metal-superconductor junc- 
tions, where Andreev reflection is suppressed at any angle 
of incidence if the Fermi wave lengths at the two sides of 
the interface are very different. The absence of reflection 
without charge conversion is a consequence of the con- 
servation of chirality (= sublattice index) by Andreev 
reflection: At normal incidence the incident electron and 
the reflected hole move on the same sublattice, while the 
reflection without charge conversion would require scat- 
tering from one sublattice to the other. The same con- 
servation of chirality is responsible for the perfect trans- 
mission of normally- incident Dirac fermions through a 
potential barrier |lj, Q . 

The differential conductance of the NS junction follows 
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FIG. 3: Differential conductance (normalized by the ballis- 
tic value go = ANe^ /h) of the interface between normal and 
superconducting graphene, for the case of small and large 
Fermi wave length \f in the normal region (relative to the 
coherence length ^ = hv/Ao in the superconductor). The 
electron-hole conversion is predominantly retro-reflection for 
Af <?C C (dashed curve), and predominantly specular reflec- 
tion for Xf 3> £, (solid curve). For eV < Ao the two curves 
are each others mirror image (when plotted versus V^). For 
eV 2> Ao both curves tend to (4 — n)go. 



from the Blonder-Tinkham-Klapwijk formula p^ . 



a/ 

dV 



1^ = 9oiV) 



7r/2 



'l-\r{eV,a)\^ + \rA{eV,a)\^) 



X cos a da, 



(13) 



9oiV) = ^iV(ey), Nie) = i^^^. (14) 

h Trnv 

The quantity go is the ballistic conductance of N trans- 
verse modes in a sheet of graphene of width W (each 
mode having a four- fold spin and valley degeneracy) . We 
assume N ^ 1, disregarding here the threshold effects 
that occur when N becomes of order unity [13, [la| ■ All 
integrals can be done analytically. The results are plot- 
ted in Fig. O for the two opposite regimes Xp ^ £, and 

Xp^C 

The differential conductance has a singularity at eV = 
Aq, as usual for an NS junction [l^. For eV ^ Aq 
we find dl/dV — > (4 — TT)go ~ 0.86 f/o, somewhat below 
the ballistic value due to the mismatch of Fermi wave 
lengths at the two sides of the interface. The subgap 
conductance, in contrast, exceeds 50 because of Andreev 
refiection. The ratio {dI/dV)/go varies between 4/3 and 
2 for both retro-reflection and specular Andreev reflec- 
tion, but the direction of the variation is inverted in the 
two cases. The difference between the solid and dashed 
curves in Fig. |3| is a unique observable signature of the 
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FIG. 4: DifTerential conductance of the NS interface for 
Sf/Ao = 0.01,0.1,0.5,1,10 (solid curves labeled b,c,d,e,f, 
respectively). The dashed curves are the asymptotes for 
Ef/ Ao ^0,00 (labeled a,g, respectively). Notice that these 
curves are plotted versus V, instead of versus V^ as in Fig.|3| 



type of Andreev reflection one is dealing with. 

In experiments it may be difficult to reach the regime 
Ef ^ Ao, so it is of importance to also consider 
the regime of comparable Ep and Aq, in which retro- 
reflection crosses over to specular Andreev reflection. 
The differential conductance in the crossover regime is 
plotted in Fig.^ (See the Appendix for the calculation.) 
It approaches the two limiting behaviors shown in Fig. 
Ofor Ep ^ Aq or Ep ^ Aq. The crossover from one 
limiting curve to the other is highly nonuniform. In the 
limit V ^ one has g^^dl/dV -^ 4/3 for any finite ra- 
tio Ep/Aq. For Ep < Ao the differential conductance 
vanishes identically at eV — Ep, because when s = Ep 
there is no Andreev reflection for any angle of incidence. 
These two conditions together imply a drop oigQ^dl/dV 
from 4/3 to as eV increases from to Ep < Ao. The 
drop becomes very rapid if Ep ^ Aq. All of this should 
be unambiguously observable. 

In conclusion, we have shown that Andreev reflection 
in graphene is fundamentally different from normal met- 
als. Close to the Dirac point (at which conduction and 
valence bands touch), an electron from the conduction 
band is converted by a superconductor into a hole from 
the valence band. The inter-band electron-hole conver- 
sion is associated with specular reflection, instead of the 
usual retro-reflection (associated with electron-hole con- 
version within the conduction band). This is but the 
flrst example of an entirely new phenomenology to ex- 
plore, regarding the interplay of superconductivity and 
relativistic quantum dynamics. We have demonstrated 
how the conductance of a single normal-superconductor 




hv\q\ > Ep 




FIG. 5: Dispersion relation in the normal state, for two 
choices of transverse wave vector: Tiv\q\ < Ef (upper panel) 
and hv\q\ > Ef (lower panel). Electron states ^"^ lie on the 
red curve, while hole states ^'' lie on the blue curve (solid 
for conduction band, dotted for valence band). The index ± 
is the sign of the slope, which equals the sign of the velocity 
in the a;-direction. 



interface (NS junction) is drastically changed by the tran- 
sition from retro-reflection to specular Andreev reflec- 
tion. We anticipate more surprises in connection with 
the Josephson effect for an SNS junction. 

Discussions with M. Titov are gratefully acknowledged. 
This research was supported by the Dutch Science Foun- 
dation NWO/FOM. 



APPENDIX A: CALCULATION OF THE 
ANDREEV REFLECTION AMPLITUDES 

1. Scattering states in N 

To calculate the Andreev reflection at the NS interface 
in graphene we first construct a basis of scattering states 
in the normal region x > 0. These are solutions of the 
DBdG equation with A = [/ = 0, 



P ■ (T Ep 





Ep 




-pa] 



(Al) 



where p ■ a = —iTiv{axdx + Uydy). 

At a given energy e and transverse wave vector q we 



have up to four basis states (see Fig.|Sl), 



* 



'exp(— ia/2)N 
,_^_ _ exp(i9y + ikx) [ exp(ia/2) 






/ exp(iQ;/2) 



^/cosa 



$<= 



eicp{iqy — ikx) 
^cosa 




^/i+ ^ cxp(zg^/ + ik'x) I 

%/cosa' exp(-ia72) 

\- exp(m72)y 

/ 

^,,_ ^ exp(^gy - zfc'o;) 

Vcosa' exp{ia'/2) 

\exp{-ia'/2)/ 



with the definitions 



a — arcsin[?iug/(e + Ep)], 
a' — arcsin[7iw(7/(e — Ep)], 
k = {hv)'^^{e + EF)cosa, 
k' = {hv)-^ {e - Ef) cos a'. 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 

(A7) 
(A8) 
(A9) 



The angle a G (— 7r/2,7r/2) is the angle of incidence 
of the electron (having longitudinal wave vector fc), and 
a' is the reflection angle of the hole (having longitudi- 
nal wave vector fc'). The signs are defined such that 
for retro-reflection a' , k' both have the opposite sign as 
a, k (this happens for e < Ep), while for specular reflec- 
tion a', k' and a, k have the same sign (this happens for 
e > Ep)- With this sign convention the states vJ/'^+j vI/'*+ 
move in the +x direction (away from the NS interface), 
while vj/*^^, \|/'*^ move in the —x direction (towards the 
NS interface). The factors l/-\/cosa and l/vcoso/ in 
Eqs. (|A2l - IA5p ensure that all four states carry the same 
particle current. 

There is no Andreev reflection beyond a critical angle 
of incidence a^, defined by 



E, 



For I a I > ac one should take 



Ep 



sign (a) ( — sign (e — Ep) — i arcosh 



2. Scattering states in S 



(AlO) 



sma 



smac 



(All) 



In the superconducting region a; < we need solutions 
of the DBdG equation 



p ■ a- - Uo - Ep Aoe*"^ 

Aoe-^"^ Ua + Ep-p-a] 



V J 

(A12) 




FIG. 6: Dispersion relation in the superconducting state, 
according to Eq. ||HJ. For e < Ao there are no propagating 
waves in the superconductor. 



that either decay as x —^ — cxd (for e < Aq) or which prop- 
agate in the —x direction (for e > Aq). The dispersion 
relation (jSJ is shown in Fig. |S1 

In the relevant regime Uq + Ep ;:^ Aq, e these solutions 
take the form 

(exp(-i^) \ 
'^expKf^ ,(A13) 
exp(i7 — iip)/ 

(exp(i/3) 
"'expH^V^^ 
— exp(— Z7 — iq 

(A14) 
The parameters /3, 7, fco, k are defined by 



f3 



arccos(e/Ao) if e < Aq, 
—z arcosh (e/Ao) if £ > Aq, 

7 — arcsin[?iwg/(t/o + Ep)], 

fco 



K = 



v/(C/o + Ep)yihv)^ - g2 
{Uo + Ep)An 



(M'fco 



■ sin/3. 



(A15) 

(A16) 
(A17) 

(A18) 



For e ^ Aq the state \p'^+ represents a hole and vj/'S'- 
an electron, both propagating in the —x direction. For 
e < Aq these states are coherent superpositions of elec- 
tron and hole excitations in the superconductor, expo- 
nentially decaying as x —^ —00. 

In the main text we assumed that Uq ^ Ep,e, mean- 
ing that the Fermi wave length X'p in the superconduct- 
ing region is much smaller than the wave length Xp in 
the normal region. Since \q\ < Ep/fiv, this regime of a 
heavily doped superconductor corresponds to the limits 
7 ^- 0, fco ^ Uo/hv, K -^ {Ao/fiv)sm(3. The states 



T3|) and HA14p thus simplify to 
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s+ 



* 



exp(— i/3)N 

exp( — Z(/ 

exp{—iq 

exp(i/?) 

— exp(z/3) 
exp(— 10) 

- exp(— i(/)) 



S— iqy — ikQX-\-KX 



(A19) 



(A20) 



3. Reflection matrix 

To complete the construction of the scattering state we 
need to match the states in N and S at the NS boundary, 
by demanding continuity at a; = 0. The reflection am- 
plitudes for an incident electron are obtained by solving 



J 



*" 



r\S<^+ + ta* = a*" 



6* 



s- 



(A21) 



at x = for r,rji,a,b. The reflection amplitudes for an 
incident hole are likewise obtained by solving 



^''- +/*'*++ r^i^'^^ 



'lT(S'+ ^ ^'lT(S'- 



a* 



&'*■= 



(A22) 



The four reflection amplitudes r, r' ,rA, r'j^ together form 
the reflection matrix 



.VA r 



(A23) 



of the NS interface. 

We give the results in the regime A^ ^ Af of a heav- 
ily doped superconductor [taking Eqs. (jA19IA20|) for the 
states in S]: 



TA 



^A 

r' 
X 



e *"*X ^Vcosacosa' if |q;| < ac, 
if |a| > ttc, 



r = X"i(-cos/3sin[(a'-|-a)/2)]-f isin/3sin[(a'-a)/2]), 



e'^^A, 



X"^ (cos /3 sin[(a' + a)/2)] + i sin/3 sin[(a' - a)/2]) , 
cos (3 cos[(a' — a)/2)] + z sin /? cos[(a' + a)/2] . 

I 



(A24) 

(A25) 
(A26) 
(A27) 
(A28) 



One can verify the unitarity condition {TZTZ')nm — Snm 
for e < Aq. (For £ > Ag the reflection matrix TZ is not 
unitary because there is also transmission of electrons 
and hole excitations into the superconductor.) One can 
also check that the two limits Ep ^ e (when a' = —a) 



and Ep <^ e (when a' = a) of Eqs. I|A24|I and l|A25|l re- 
duce to Eqs. l| ^12|l in the main text. The full expressions 
given above were used to calculate the crossover curves 
plotted in Fig. 2| 
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